pacman::p_load(ggtree,dplyr,stringi)
setwd("E:/projects/R_draw_sc/test/hca/data")
data <- read.delim("hca.list",row.names = 1,check.names = FALSE)

col <- c("z13-15","z14-8","yh2")
row <- c("protein-5","POD-10","SOD-5","SOD-10")

data <- data[,match(col,colnames(data))]
data <- data[match(row,rownames(data)),]
data %>% t() %>% dist() %>% hclust()%>% plot()
